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Abstract 

The  present  paper  deals  with  the  one -dimensional 
propagation  of  disturbances  in  an  invlscid  conducting 
fluid  of  finite  magnetic  Reynolds  number.   The  following 
conspicuous  property  is  shown:   The  basic  equations  are 
not  hyperbolic  but  nevertheless  the  slow  wave  has  a  domain 
of  dependence  determined  by  the  sound  velocity. 
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The  Effect  of  Finite  Conductivity  on  the  Propagation  of 
Hydromagnetic  Slow  Waves 

1 .   Introduction 

The  effect  of  finite  conductivity  on  the  propagation  of 
hydromagnetic  waves  of  small  amplitudes  [l]  has  been  investi- 
gated by  Cole  [2]  and  Whltham  [3]  for  the  transverse  case  in 
which  the  magnetic  field  is  perpendicular  to  the  direction  of 
wave  propagation.   The  present  paper  deals  with  one  dimensional 
propagation  in  a  general  configuration  where  the  magnetic 
field  is  oblique  to  the  direction  of  wave  motion.   Even  if 
the  fluid  is  Invlscid,  the  finite  conductivity  makes  the  basic 
equations  parabolic  so  that  in  the  usual  sense  the  wave  front 
proceeds  with  infinite  speed  and  the  existence  of  the  domain 
of  dependence  will  be  lost  completely.   In  the  present  case 
this  is  true  for  the  fast  wave;  the  slow  wave,  however,  has  a 
domain  of  dependence  determined  by  the  sound  speed,  the  wave 
front  of  which  proceeds  with  finite  speed.   In  Section  2  the 
dispersion  relation  is  derived,  in  Sections  3  and  4  the  above 
mentioned  peculiar  property  of  the  slow  wave  is  shown  on  the 
basis  of  the  Laplace  transform,  and  in  the  last  section  the 
effect  of  the  Hall  current  is  taken  into  account. 

The  author  would  like  to  express  his  gratitude  to  Dr. 
C.  S.  Gardner  for  his  advice. 
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2 .   The  Basic  Equations  and  the  Ass-umptions . 

In  this  paper  we  consider  the  one  dimensional  propagation 
of  disturbances  In  a  quiescent  conducting  fluid  In  a  constant 
external  magnetic  field  H   so  that  all  the  quantities  will  be 
assumed  to  be  the  functions  of  x  and  t   only.   The  effect 
of  viscosity  will  be  neglected  throughout,  while  that  of  finite 
conductivity  will  be  taken  Into  account;  however.  In  small 
disturbances  the  joule  heating  resulting  from  the  finite  con- 
ductivity can  also  be  neglected  as  higher  order,  hence  we  may 
assume  that  the  flow  Is  Isentroplc. 

The  coordinate  system  Is  so  Introduced  that  the  z -component 

of  H  ,  H   ,   Is  zero.   The  disturbances  of  the  fluid  velocity, 
o   zo  ^ 

the  density,  and  the  magnetic  field  are  denoted  by  v  ,  v  ,  v  , 
p,  and   H  ,H    respectively,  while  the  unperturbed  constant 
density  and  the  corresponding  sound  velocity  are  denoted  by  p 
and  a   respectively.   By  virtue  of  the  dlvergenceless  condi- 
tion of  the  magnetic  field,  the  disturbance   H,   can  be  put 
equal  to  zero.   Neglecting  the  terms  quadratic  with  respect  to 
the  disturbed  quantities  we  have  the  linearized  system  of 
equations 


(2.1) 
(2.2) 


•X                                              SV 

3t            "PQ     ax 

PQ     ^t                "     ^0    ^    ■ 

HTT      yo   dx 

In  this  paper  we  employ  Gaussian  units  and  m-  denotes  the 
magnetic  permeability,  c  the  velocity  of  light  In  the  vacuum. 
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X 


(2-3)  Po3t^=feH^„5 


dx 
SH  av  2      d^H 

(2.5b)  -vT^    =     H  ^r-^    +     -r-^ ^ 

^      -^    ^  dt  xo   dx  4Tri-La    -.    2 


in  which  the  ordinary  Ohmic  law  for  the  electric  field  vector 
E  and  the  current   j , 


(2.6)      a(E  +  ^  [v  XH])  =  J 


is  assumed.   Introducing  the  dimensionless  quantities  denoted 
by  a  prime  through  the  equations 


x  =   Lx',         t    =   TtS  H=H^H',  p=pp', 

y         yo  y  "^  o*^ 


r  XI 

"•o        ■'^"'"      '       I^Trp .      ~    'xo        'xo 


a      =   LT    ^,       V^^^     5   b.„    =   b,'.„    LT"^ 


) 


^           yo          yo  ^T^ 

and   eliminating     v        and  v        from  Eqs .     (2.1)    -    (2.4),    we   get 

^  y 

.,2              .2    \  o    ^^H 

dt'2     ax'2  y  y°  ax '2 


<'•'"     Ut.2       °    d.'^l   y     ax'2      "■   ax'^at' 


in  which  Rm  stands  for  the  magnetic  Reynolds  number,  given 
a  L  ^Tr[io 

by   Rm  =  p 

c 
In  what  follows  all  the  primes  of  these  dimensionless 

quantities  will  be  omitted  unless  otherwise  stated.   The 

system  of  equations  (2.7  a,b)  shows  the  coupling  between  the 

density  compression  and  the  magnetic  disturbance,  and  in  the 

limit  as   Rm  — >  oo,  leads  to  the  fast  and  slow  waves  [l],  while 

equations  (2.5)  represent   the  transverse  wave  [l]  and  are 

still  separated  from  the  remaining  equations  even  for  finite  Rm. 

In  order  to  obtain  the  dispersion  relation  resulting  from 

equation  (2.7)  we  put 

(P    \^eP^+^^ 

^  y ' 

Inserting  the  above  expression  into  equations  (2.7  a,b)  and 

2     2         2 
denoting  b  ^,  b,^^  and  b^  by  r^,  r  and   r  respectively, 
xo    yo        o       X    y 

we  obtain  the  secular  equation 


2 
(2.8)      q^  =  S-  (^  +  Rm"^p)~^[l  +  r  +  Rm"^p 


X 


±^(1  +  r  +  Rm  ^p)^-  ^(r  +  Rm'^p)  ]   . 

We  now  introduce  the  complex  phase  velocities   A   and  A^ 
through  the  equations 

(2.9)     >^^^2  =  ^1  ^^x+  Rm'^p)"^h  +  r  +  Rm"^p 

+  /(I  +  r  +  Rm'^p)^-  4(r^+  Rm"^p)  j]^/^ 


In  which  subscripts  1  and  2  correspond  to  the  minus  and  the 
plus  signs  in  the  right  hand  side  respectively;  it  is  easy 
to  see  that  A  and  A   tend  to  the  fast  and  the  slow  speed  [l] 
as  Rm  approaches  Infinity.   Then  for  any  complex  value  of  p, 
q   is  given  by  the  equation. 


(2.8).  q.  =±^,       (i=i,2) 


3-   A  Boundary  and  Initial  Value  Problem 

We  now  consider  an  initial  value  and  boundary  value 
problem,  defined  as  follows.   The  boundary  condition: 
at   X  =  0,   U  is  given  as  a  function  of  time  in  the  following 
way, 

(A.l)  f  =   ]j    (t)  for   t  >  0 

U(o,t)  < 
(A. 2)  =  0        for   t  <  0  ; 


the  initial  condition: 

at   t  =  0,   U  and  U,   are  zero  for  positive   x,   namely, 

U(x,o)  =  0 
(B)  \    for  X  >  0 

U^(x,o)  =  0 

where  U   is  a  column  vector  defined  by 
(3.1)     U  = 


{;) 


The  solution  for  the  positive  value  of  x,   satisfying  the 
above  conditions  (a)  and  (B)  can  be  given  by  the  Laplace 
transform, 

(3.2)      U(x,t)'-  2?r  f  ^^^^'^'^   U(p)dp 

In  which  q   Is  given  by  q.   of  equation  (2.8).   (The  com- 
plete solution  Is,  of  course,  the  superposition  of  the  solutions 
corresponding  to  the  various   q. 's.)   The  contour  L  Is  the 
Bromwlch  path,  with  Re  p  >  0,   to  the  right  of  all  singularities 
of  the  Integrand.    Deforming  the  Integration  path  to  encircle 
the  positive  half  plane,  one  can  see  Immediately  that  the 
boundary  condition  (A. 2)  Is  automatically  satisfied,  while  the 
Initial  condition  (b)  requires  that   Re  q  must  be  negative 
for  Re  p  — >  oo.   Since  for   |p|  — >  «>,   A.  's  behave  as 

(3.3a)         ^1^^^  (^"i~^  P)~^  ^ 
(3.3b)         A^^  ->  1 

the  above  requirement  admits  among  the  q. 's  of  equation  (2.8) 
those  with  negative  sign, 

(3.4)  q.  =  _  ^   ,     (i  =  1,2)  . 

1 

The  solution  Is  the  superposition  of  the  elementary  solutions 
with  q. 's  given  by  equation  (3.4)  which  clearly  expresses  the 


These  singularities  comprise  not  only  those  of  U(p)  but  also 
the  branch  points  of  q  which  are  given  by 

Rm~  p  =  -r^,     Rm~"'"p  =  1  -  r  t  21 V?"   . 
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outgoing  wave.   Prom  equation  (2.7a)  It  follows  that  the 
solution  takes  the  form, 

(3.5)  U(x,t)  =  2?i  _[  S(p,x)S-^  U^(p)eP^  dp 

in  which  S(p,x)  is  defined  by 


Ye  7   6 

y  .         y 


(3.6a)         S(p,x)  - 


S    is  equal  to  S   (p,o)   and  has  the  representation 


(3.6b) 


/ 

/- 

-I 
-I 

1 

-r'^- 
-.1^^^ 

1 

0 

^2         ^2 
^1    -    ^2 

1 

T.2         -.2 

while  U  (p)   is  the  Laplace  transform  of  U  (t), 

00 
(3.6c)         U^(p)  =  f  e-P^  U^(t)dt. 

Equation  (3«5)  can  also  be  written  as 

-,  f  pt+q-,x       pt+q„x 

(3.7)  U(^'^^=2^J   (^1^    ^  +  A2  e    2  )^p 


if  the  column  vectors  A   and  Ap   are  introduced  through 


the  equations 
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(3.8a)    A^  =  (S-^  U^(p))l 


X^I 


(3.8b)     A„  =  (S"-^  U  (p)) 


2    ^  o   "o^^^^2 


/y 


\  ^2  -  1 


where   (S   U  (p)).,  (i  =  1,  2)  denotes  the   i    component 
of  the  column  vector   S   U  (p). 
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4 .   The  Domain,  of  Dependence  of  the  Slow  Wave 

Let  us  first  evaluate  the  integral  for  the  slow  mode 
pt+qgX 
characterized  by  e        in  equation  (5-7). 

From  equation  (j.Jb)  it  follows  that   q^   tends  to  -p, 

consequently   pt+q„x,   to   p(t-x),   as   Re  p   approaches  +oo. 

Hence  deforming  the  path  L  to  encircle  the  positive  half 

2   (     Pt  +  qg^ 
plane,  one  can  see  easily  that  the  integral  ^^  1   ^p^        '^P 

vanishes  if   t-x  <  0.   This  shows  the  striking  feature  that 
even  for  the  finite  value  of  conductivity  the  slow  disturbance 
has  a  domain  of  dependence  and  does  not  propagate  with  a  velo- 
city larger  than  the  sonic  one.   On  the  other  hand,  the  fast 
mode  behaves  in  an  usual  sense;  namely,  from  equation  (3.5a) 
q-,   behaves  as   -/p  as  Re  p  tends  to  infinity  and  the  same 
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deformation  of  path  does  not  lead  to  the  existence  of  a  domain 

of  dependence. 

It  seems  worth  while  to  check  these  conclusions  by  the 

method  of  steepest  descent,  for  the  asymptotic  behavior  of 

large   t.   In  what  follows  we  consider  the  case  that  the  external 

magnetic  field  Is  neither  parallel  nor  normal  to  the   x  axis 

so  that  both   r   and   r   differ  appreciably  from  zero.   The 
X         y 

special  configuration  r  =  1  will  also  be  excluded  from  our 
consideration.   Let  us  for  a  while  assume  that   Rm  p   is  small; 
then  we  may  expand   q.  =  -  p/-^.   as  a  power  series  in  terms  of 
Rm  p, 

(4.1)      q.  =  -  ^  [1  +  ^  E  +  0{e^)]    , 

where   e   stands  for  Rm  p  and  only  the  zeroth  and  the  first 
order  terms  in   e   are  written;  A.  and  A^   are  the  fast  and 
the  slow  speed  given  by  [l]. 


(4.2)      Ai  =  ^1  (Rm-1  =  ^)    jjZ^V^^f^pLHilL 


2  2 


and  a.   is  defined  by 
1 


(4-3)     a^  =  -(Ai)'^  (1  t  ^"^       )  • 

2       2         /(l-r)^  +  4r 

(The  subscripts  1  and  2  correspond  to  the  +  and  -  signs  in  the 
right  member,  respectively.)   It  should  be  noted  that  the  a.'s 
are  always  negative  and  finite.   In  deriving  equation  (4.1)  we 
assumed  not  only   e  <  1   but  also 
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(4.3) '    ea^  <  1  . 

Using  the  expression  (4.1)  for  the  Integral  (3-7)   we  have  the 
saddle  point  given  by  the  equation 

^^1  X 
dp  t 


namely. 


-1        Ai  I  -  1 


(4.4)     Rm  "p.  =    ^ 

1 

If  x/t ,  kept  constant.  Is  sufficiently  close  to  .\..,   not 
only  Rm  p.   but  also  Rm  p. a.  can  be  small  so  that  the 
expansion  used  above  Is  justified.   Therefore,  one  can  conclude 
that  the  main  disturbance  of  the  fast  and  the  slow  mode  Is 
given  by  the  fast  and  the  slow  speeds  /i-,   and  /- „ .   This 
result  Is  valid  for  the  fast  wave  In  any  case;  however,  for 
the  slow  wave  some  care  must  be  taken.   Since  there  Is  an 
Inequality  [l] 


(^•5)      Ag  <  mln  (1,  b^), 

we  may  conclude  that  the  main  disturbance  of  the  slow  mode 

propagates  with  the  speed  ,  ip   and  decreases  tc.ard  the  boundary 

determined  by  the  sound  speed,  provided  that  the  Alfven  velocity 

Is  smaller  than  the  sound  velocity,   bo  <  1;  however  If  bo  >  1 

we  must  be  more  careful.   For  Instance,  let  us  suppose  In  this 

case  that   r   approaches   zero  so  that  /i„   tends  to   1,  then 
y  2  ' 

Qp   also  approaches  zero  since   r  >  1;  consequently  the 


13  - 


assumption  (^.3) 'ceases  to  be  valid.   From  the  previous  result 
one  could  expect  that  an  adequate  method  would  give  the  dis- 
turbance rapidly  decreasing  toward  the  wave  front   x-t  =  0, 
vanishing  there  exactly.   If   r   is  zero,  we  can  see  easily  that 
equations  (2.1) -(2. 5)  become  completely  separated  into  the  sonic 
and  the  Alfven  equations  and  the  slow  wave  becomes  the  sound 
wave . 


5.   The  Effect  of  the  Hall  Current 

In  a  plasma,  the  Ohmic  law  must  be  modified  as  follows  [^] 


(5.1)*    a(E  +  ^VXH  +  —  VP)  =  J+  ^  j  >^  H 
\^      I  \  2         nee        mc"^ 

e 

where   e   is  the  electric  charge,   m   the  electron  mass,   t 
the  collision  frequency,  and  n  the.  density  of  plasma. 

Assuming  p  =  p/2 ,   and  using  equation  (5.I)  for  the 
linearized  one  dimensional  propagation  we  have 

^2     -^2       ^   B^H 

(5-2)    (f;?  -  r^)p  -  ''yo  T^  =  °  ■  ■ 

dt     dx        ^   dx 

x2      -,2       ^2^       ,  B^H  T  a^H 

(5.3)      (^  ^^)H  -  ^  -  Rm"^  —^  -    (C  tan  ^)Rm~-^ 


St     b^ax     ^    ax  Sx  at  Sx  at 


=    0 


We  neglect  the  inertia  of  electron  (the  plasma  oscillation) 
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(5.4)     (^  -  b2^  ^)H^-  Rm-1  -^  +  (C  cot  ^)Rm-l  -^  =.  0 
5t        Sx   '^       Sx  St  dx  St 


in  which  all  the  quantities  x,  t,    p,  H   and  H   denote  the 
dimensionless  ones  used  so  far,  and  i/   is  given  by  the  equation 


while  C  is  defined  by 


tan  tA  =  H  /H 

zo   yo 


M-eH 
^      xo 

^    m  c 
e 


If  the  mean  free  path  of  the  electron  is  sufficiently 
smaller  than  the  Larmor  radius  of  the  electron,  then  C  is  small 
and  the  effects  of  the  Hall  current  represented  by  the  last 
terms  in  equations  (5-2)  and  (5-3)  can  be  neglected,  so  that 
the  magneto-accoustic  wave  and  the  transverse  wave  are  decoupled. 
However,  even  if  C,   is  not  negligibly  small  it  can  be  shown  that 
there  still  exists  a  mode  of  propagation  with  the  domain  of 
dependence  given  by  the  sound  velocity. 

The  secular  equation  corresponding  to  equation  (2.o)  becomes 

(5.5)     [(r^  +  ef   +    {Cef]q^-   p^[(r^+  £)(2  +  r  +  e)  +  U^fl^i"^ 

+   p  (1  +  r+  r^  +  2e)q^-  p  =  0. 
In  the  limit  |p|  — >•  oo,  the  above  equation  may  be  approximated  by 
(C^+  l)£2(q/p)6_  (^2^  i)e2(q/p)^  +  2e{q/pf    -1  =  0 

hence  we  still  have  one  root  q  which  behaves  as  q  --'  -p  for 
Re  p  — ^  00,  consequently  the  mode  corresponding  to  this  root 
has  a  domain  of  dependence  associated  with  the  sound  velocity. 
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